Carbon nanotubes in electric and magnetic fields 
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We derive an effective low-energy theory for metallic (armchair and non-armchair) single-wall 
nanotubes in the presence of an electric field perpendicular to the nanotube axis, and in the presence 
of magnetic fields, taking into account spin-orbit interactions and screening effects on the basis of 
a microscopic tight binding model. The interplay between electric field and spin-orbit interaction 
allows us to tune armchair nanotubes into a helical conductor in both Dirac valleys. Metallic 
non-armchair nanotubes are gapped by the surface curvature, yet helical conduction modes can 
be restored in one of the valleys by a magnetic field along the nanotube axis. Furthermore, we 
discuss electric dipole spin resonance in carbon nanotubes, and find that the Rabi frequency shows 
a pronounced dependence on the momentum along the nanotube. 

PACS numbers: 73.63.Fg, 72.25.-b, 75.70.Tj, 85.75.-d 



I. INTRODUCTION 

The last two decades have seen remarkable progress 
in the experimental techniques to fabricate and analyze 
high purity carbon nanotubes (CNTs)ii~— This progress 
has paved the way for using CNTs for electron and, in 
particular, electron spin based effects that are of interest 
for quantum information processing and spintronics£~— 
For the latter, the spin-orbit interaction (SOI) plays a sig- 
nificant role as it allows the spin manipulation by electric 
fields. The purity of CNTs has by now advanced so far 
that indeed SOI effects can be observed^— 

In this paper, we investigate the SOI in metallic single- 
wall CNTs in the presence of external electric and mag- 
netic fields. We provide an extensive discussion of electric 
field screening in CNTs and of resonant spin transitions 
via ac electric fields, namely, electric dipole spin reso- 
nance (EDSR). Our starting point is a tight binding de- 
scription including all second shell orbitals of the carbon 
atoms, from which we derive an effective low-energy band 
theory. 

In the metallic regime, SOI effects by an external elec- 
tric field can only be expected if the field is applied per- 
pendicularly to the CNT axis. Without the electric field, 
the SOI in a CNT has been studied before^— and is the 
result of the orbital mixing caused by the curvature of the 
graphene sheet wrapped up into a cylinder. This leads 
to different geometric conditions than for flat graphene, 
which has accordingly a different SO L 30 : 31 It is also very 
different from the SOI found in typical one-dimensional 
semiconductor wires because of the rotational symmetry 
of the CNT. In a semiconductor wire, the usual Rashba- 
SOI with a well-defined spin-precession axis arises from 
the specific asymmetric electric environment caused by 
the confining potentials. Such an axis is absent in the ro- 
tationally invariant CNT unless it is reintroduced by the 
application of an external transverse electric field. As 
we discuss in detail below, this leads indeed to a Rashba- 
like SOI. However, in contrast to the semiconductors, the 
CNTs have a hexagonal lattice structure with two carbon 
atoms per unit cell. The resulting band structure has a 



nearly vanishing density of states at the charge neutral- 
ity point (Dirac point), and roughly a linearly increasing 
density of states at low energies away from this point. 
This means, first, that the high energy states affect the 
low energy properties in a more pronounced way than in 
a semiconductor. Second, the screening of the external 
electric field becomes less effective, and the distribution 
of charges on the tube surface can become complicated. 
As the screened field affects the low-energy physics as 
well, it requires specific investigation. It turns out here 
that for a quantitative understanding we need to start 
indeed from the full band structure based on the lattice 
description, and not from the effective low-energy Dirac 
theory. 

The combination of these effects allows us to under- 
stand how precisely SOI and especially the SOI parts 
induced by the external electric field affect the system 
properties. The larger number of external and internal 
degrees of freedom as compared to semiconductor wires, 
such as field strength, sublatticc or Dirac valley index, 
CNT radius, chirality, and chemical potential, allows for 
an extensive tunability of the SOI-induced system prop- 
erties. For instance, in armchair CNTs it allows us to 
obtain in both Dirac valleys helical, spin-filtered conduc- 
tion modes, in which opposite spins are transported in 
opposite directions. Similar helical states occur in SOI- 
split quantum wires^r— and at the edges of topologi- 
cal insulators^ and can be used, for instance, as spin 
filter o 39 ' 41 or as Cooper-pair splitters42 In addition, if 
such a helical conductor is brought into proximity to 
a superconductor, it allows for the realization of Majo- 
rana end states. This has attracted much attention very 
recently^ 3 -"— mainly because these Majorana states may 
be used as fundamental quantum states for topological 
quantum computation^ 

It should be noted that in semiconductor quantum 
wires the helical modes are realized only under condi- 
tions of an external magnetic field, with the consequence 
that time-reversal symmetry is broken and the spins at 
opposite conduction band branches are not truly antipar- 
allel. In contrast, in armchair CNTs the helical modes 
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can be obtained in an all-electric setup, and they are 
perfectly polarized in the sense that only the spin ex- 
pectation value S y in the direction perpendicular to the 
electric field and the CNT axis is nonzero, even though 
we generally have \(S V ) \ < 1 within the helical bands. 

For metallic but non-armchair CNTs, two more terms 
appear in the low-energy Hamiltonian. First, an orbital 
curvature induced term, which opens gaps at the remain- 
ing helical zero-energy modes. Second, a SOI term, which 
plays the role of an effective Zeeman field transverse to 
the S y polarization, yet with opposite sign in each Dirac 
valley. While these two terms primarily destroy the heli- 
cal conduction modes, their presence can be turned to an 
advantage: by applying a magnetic field along the CNT 
axis, the effective Zeeman field in one of the valleys can 
be suppressed, while it doubles the gap in the other val- 
ley. Hence the helical conduction modes are restored to 
high quality in one of the valleys, while conduction in the 
other valley is suppressed completely. 

To conclude, we provide a microscopic description of 
electric dipole spin resonance^ - — (EDSR) in CNTs. 
An ac electric field perpendicular to the CNT axis cou- 
ples to the electronic spin via the SOI. However, this cou- 
pling comprises an additional sublattice coupling that is 
absent in semiconductor setups. It causes a significant 
momentum dependence of the resonant Rabi frequency 
of the EDSR experiment. The further application of a 
static electric field, perpendicular to the CNT axis and 
to the ac field, lifts the spin degeneracy of the bands in 
an way analogous to a static Zeeman field. This allows 
us to propose an all electric realization of Rabi resonance 
experiments, in which the electric fields replace the static 
and timc-dcpcndcnt magnetic fields. 

The paper is structured as follows. In Sec. [II] we in- 
troduce the tight-binding model of the hexagonal carbon 
lattice including SOI and the effect of electric and mag- 
netic fields. The 7r and a bands are discussed in Sec. 
IIII1 and we identify the important ira band hybridization 
couplings that have an impact on the low-energy physics. 
This allows us to derive the effective low-energy theory 
for the 7r bands in Sec. IIV1 The partial screening of the 
external electric field is investigated in Sec. [Vj where we 
discuss in particular the validity of linear response and 
the influence of the er bands. The implications of the ef- 
fective low-energy theory are analyzed in Sec. I VII with 
special focus on helical modes, valley suppression and the 
role of further external magnetic fields. Turning to the 
dynamical response of the CNTs, in Sec. IVII| we investi- 
gate resonant spin transitions by an ac electric field. The 
final section [Villi contains our conclusions. 



II. MODEL 

Our calculations are based on a tight-binding model 
describing the second shell orbitals of the carbon atoms in 
CNTs.— It takes into account the spin-orbit interaction 
and external electric and magnetic fields. The Hamilto- 



nian we start from consists of four terms 

H = i?hop + Hso + He + Hb, (1) 

where i/hop describes the hopping between neighboring 
tight-binding orbitals, Hso describes the on-site spin- 
orbit interaction, and He and Hb describe external elec- 
tric and magnetic fields, respectively The Hamiltonian 
is expressed in second quantization by the electron an- 
nihilation operators c n ^ iAii A, where the two-dimensional 
integer vector n = (711,712) labels the unit cells of a 
honeycomb lattice, ( = ±1 labels the sublattice A/B, 
and A = ±1 labels the spin, /j = s,p ri p tl p z denotes 
the second shell orbitals with fi = s the s orbital and 
H = p ri pt,p z the p orbitals. The p r orbital is perpendic- 
ular to the surface of the CNT and essentially constitutes 
the 7r band. The pt. z orbitals are tangential (see Figs. Q] 
and[2]) and, together with the s orbital, build the a bands. 
These orbitals span a local coordinate system. 

If the cylinder surface on which the carbon atoms are 
placed in a CNT is projected onto a plane, the graphenc 
lattice with lattice vectors ai = a(t cos 9 + z sin 9) and 
a2 = a(tcos(7r/3 — 6) — zsin(7r/3 — 9)) is recovered, 
where 9 is the chiral angle of the nanotube, defined as 
the angle between the first unit vector ai and the chi- 
ral vector C (see Fig. [I]). In the following, we will 
also use the collective site index i = (711,712, C) to la- 
bel the atom positions. The position of atom i is then 
R; = niai + ii2a2 + (1 — C)L 1 /2. The three vectors 

L 1 = -^=(— tsinfl + zcosfl), (2) 
v3 

L 2 = -^=(tcos(7r/6-0) -zsin(7r/6-0)), (3) 
v3 

L 3 = -^=(-tcos(7r/6 + 0) -zsin(7r/6 + 6»)), (4) 
v3 

connect an A site with its three nearest neighbors on 
the B sublattice. These three B sites will be denoted 
by -Bi.2,3 in the following. Note that we have chosen to 
define the model directly in a coordinate system which 
will be convenient for the subsequent calculations. 

In the following, we describe each term in Eq. ([T]) and 
its derivation in detail. 



A. Hopping Hamiltonian 

The term -f/hop describes the hopping between orbitals 
on neighboring carbon sites as well as the on-site orbital 
energies. We only take into account nearest neighbor 
hopping and assume the orbitals on neighboring carbon 
atoms to be orthogonal. The hopping Hamiltonian has 
the form 

#hop = *I«V C U C J>' A + £s ^ 4a c ^a, (5) 

where runs over nearest neighbor sites, and e s = 

—8.9 eV is the orbital energy of the carbon s orbitals 
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FIG. 1. (Color online) Definition of the honeycomb lattice 
vectors and the direction of the p z and pt orbitals. The white 
arrows in the p orbitals point into the direction where the 
wave function is positive. The ellipses indicate the unit cells 
consisting of A and B sublattice atoms. The p orbital align- 
ment is equal on both sublattices. Together with the p r or- 
bital, which points out of the plane (not shown here), the 
directions of p r ,pt,p z form a right-handed set of vectors. 



where indices {/i, v} = {r, t, z} describe p orbitals and 
n = {1,2,3} denotes the neighboring B atoms. There- 
fore, the hopping matrix elements are 



C" = V pv (Pm ' P?" ~ (P? ' in )(P?" ' r )) + 
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FIG. 2. (Color online) Cross section of a CNT. The in- 
direction is along the nanotube axis. The orientation of the 
orbitals p r ,Pt as well as the local coordinate system f,t de- 
pends on the azimuthal angle ip. The s orbital is indicated by 
the dashed circles. 



FIG. 3. (Color online) Hopping between atomic orbitals of 
neighboring carbon sites in the graphene limit, (a) - (d) show 
the orbital combinations with non-zero hopping amplitude, 
(e), (f) show hoppings which are forbidden by symmetry. 





relative to the p orbital energy, The latter is set to zero 
in this paper. The amplitude i*?, , of an electron hop- 
ping from the orbital pi! on the jth site to the orbital (i 
on the ith site is a linear combination of the four funda- 
mental hopping amplitudes^ V ss = —6.8 eV, V sp = 5.6 
cV, V pp = -3.0 cV, V° v = 5.0 cV (see Fig. © with 
coefficients depending on the relative orientation of the 
orbitals fj, and \j! 'M- For the explicit calculation of t*jL,, 
the p orbitals are decomposed into components parallel 
and perpendicular to the i — j bond. These components 
can be easily expressed in terms of scalar products (see 
Fig. [4]). If the direction of a p orbital on the A atom is 
and the direction of a p orbital on the B n atom is p^™ , 
then the projections on the unit vector 1™ = L™/| L™ | in 
the direction from atom A to atom B n are 



f-L = (P 

CT II = Pu 



p)(p*"-r) 



l")(p 



B„ 



(6) 
(7) 



(Color online) The hopping element between a) p^ 
b) p^ and p^ n can be calculated by the decomposi- 
tion of the p^(pv n 



FIG. 4. 

and s B " 



orbital into two parts: perpendicular to 



the bond p^x(pf 1) and parallel to the bond p^u(p 



B. Spin-orbit interaction 

The spin-orbit interaction is modeled by projecting the 
Thomas Hamiltonian for a spherically symmetric atomic 
potential V(r)& 



rrintr 



1 dV(r) 
r dr 



L S 



(11) 



onto the subspace spanned by the second shell orbitals of 
each individual carbon atom. L is the electron angular 
momentum operator and S is the vector of the spin Pauli 
matrices with eigenvalues ±1. 
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By symmetry, this projection can be reduced to a 
form depending on a single parameter only, defining the 
strength of the intrinsic SOI 

Aso = (Pr I (dt[d r V] - d r [d t V}) \pf). (12) 

In the literature one finds a wide spread for the spin-orbit 
coupling constant Ago (see Refs. HH, H3, HI) . Note also 
that there are different conventions for the spin operators 
used in the spin-orbit Hamiltonian. In some works, Pauli 
matrices with eigenvalues ±1 are used while in other pa- 
pers, the spin operators have eigenvalues ±i. In this pa- 
per, we always use Pauli matrices with eigenvalues ±1. 
For the spin-orbit coupling constant one finds values be- 
tween 3 meV— and 20 meVp^ In this work we follow 
Ref. and use Ago = 6 meV. The projected spin-orbit 
Hamiltonian reads 

H so = *A S o ]T ^'c^,^ . (13) 

i,X.\' 
AW) 

Here, p,rj = p ri p tl p z run over the p orbitals only. e^ vr> 
is the Levi-Civita symbol. The index v = r,t, z labels 
the spin components in the local coordinate system, i.e., 
S r = S x cos ifi + S y sin ifi, S* = S y cos ipi — S x sin ipi , with 
(fii the azimuthal angle of site i (see Fig. [2]) and S x ' y ' z 
the spin Pauli matrices (with eigenvalues ±1). The spin- 
orbit energies emerging from ifgo and from the curvature 
effects are found to be much larger than the spin-orbit 
energies due to the d-orbitals^ allowing us to neglect 
the latter. 



C. Electric fields 

In this section we introduce the Hamiltonian describ- 
ing a homogenous external electric field E perpendicular 
to the CNT axis. Due to screening effects, local fields 
with complicated spatial dependence are generated by 
the rearrangement of the electron density on the lattice 
and have to be taken into account in principle. However, 
it turns out that none of these nontrivial contributions 
affect the physics in an essentialy way. 

We start from the most general electrostatic potential 
</>tot(r) = ^ext(r) + </>ind(r), which is the sum of the po- 
tential cxt coming from the homogeneous external field 
and the generally complicated potential </>i n d from screen- 
ing effects. The matrix elements (/z, i\ 4>(r) |z/, j) between 
the orbitals of two different carbon atoms i ^ j are much 
smaller than the typical hopping elements (~ eV) so that 
we can safely neglect them and consider only on-site ef- 
fects of the electric potential. 

The averaged potential ^ MM (R,) = (/x, i\ ^tot(r) i) 
for the fi orbital on site i changes the electrostatic energy 
of an electron in this orbital. Note that R^ is the position 
of the ith carbon atom in three-dimensional space. In the 
following, we assume that the dependence of WI (R^) on 



the orbital /z is negligible, i.e., (/> to t(Ri) — wl (Ri). Fur- 
thermore, we show in Sec. |V]that even with the induced 
potential 0i n d included, </>tot essentially depends on the 
azimuthal angle (pi similar to the case of a homogeneous 
field. Therefore, the diagonal matrix elements of the elec- 
trostatic potential are </>(Ri) ~ <fitot(<Pi) oc cosifi and give 
rise to the Hamiltonian 

H [ E ] = E ^tot^H^A. (14) 
z,/i, A 

For sufficiently small CNT radii R and electric fields the 
total potential on the surface of the CNT is well approx- 
imated by </>tot(<y5j) — eE*Rcosipi, where E* < |E| is the 
screened electric field and e is the electron charge. 

We note that, because 4> varies on the scale of the spa- 
tial extent of the orbital wave function and breaks the lat- 
tice symmetries, in general, the matrix elements between 
the different orbitals 4>^^ (R^) = (zt, i\ </>(r) \ fi' , i) ^ and 
transitions between orthogonal orbitals /i, p! on the same 
carbon atom are generated. We call this the tran- 
sitions in the following. The potential of a homogeneous 
external field alone gives rise to s-p r and s-pt transitions. 
But the complicated additional induced potential </>; n d ( r ) 
gives also rise to coupling between other orbitals. For this 
work, the s-p r transition is most important because of 
two reasons: 1) It is the only transition directly coupling 
7r and a bands, thus giving rise to a first order effect 
in the s-p r coupling strength. 2) Its strength is deter- 
mined by the unscreened field E and not by the screened 
E* < E. Indeed, the induced potential 0i n d drops out, 
i.e., (p r \ </>md(r) \ s) ~ 0, as <pi n d is approximately an even 
function in the radial coordinate r about the tube radius 
r = R. The Hamiltonian describing the s-p r transition is 

= -eE£ cos(ifi)c\ prX c ls x + H.c, (15) 

with the strength of the transitions being characterized 
by the integral 

£o = -J d 3 r^ s (r)z^(r) = (16) 

where ip2s and ipip^ are the hydrogenic wave functions of 
the second shell atomic orbitals, as is the Bohr radius, 
and Z is the effective nuclear charge, which for the second 
shell in carbon is Z ~ 3.2. From this we obtain £o ~ 
0.5 A. Note that our value for £ is a rather conservative 
estimate. It is roughly four times smaller than what has 
been assumed in Ref. l3ll . 

D. Magnetic fields 

The Zeeman Hamiltonian describing the interaction of 
the electron spin with the magnetic field is in the tight- 
binding model written as 

H mag = E ^BB v c\^ x S v xx ,c lfl x', (17) 
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where the index jj, runs over all orbitals, i labels the po- 
sition of the atom, and A, A' denote the spin. B v with 
v = x,y,z are the components of the magnetic field in 
the global coordinate system, which we choose such that 
the electric field is always in x direction and the CNT 
axis is along z. 

The orbital effect of the magnetic field is usually ex- 
pressed in terms of Peierls phases, multiplying the hop- 
ping amplitudes in the tight-binding Hamiltonian. For 
a magnetic field along the CNT axis, which will be of 
most interest in this paper, this will lead to an Aharonov- 
Bohm shift of the circumferential wave vector, which can 
be incorporated easily into the effective model, derived 
subsequently. 



III. 7T AND a BANDS 

In general, all four orbitals of the second shell in CNTs 
are hybridized by the hopping Hamiltonian -ffhop- How- 
ever, in the limit R — > oo, the hopping amplitudes be- 
tween the p r orbital and any of the orbitals s,pt,p z be- 
come zero, and the hopping Hamiltonian becomes equal 
to the Hamiltonian of flat graphene. The p r orbitals then 
form the 7r band and the remaining three orbitals, which 
are still strongly hybridized, form the a bands. For CNTs 
with finite _R, the ir band hybridizes with the er bands due 
to the CNT curvature, which leads to non-zero hopping 
amplitudes tfj Pt , tfj p *, tfj s . However, the ira hybridiza- 
tion is still small for realistic CNT sizes, so that it is 
convenient to partition the total Hamiltonian [Eq. (JTJ) ] 
as 

H = H n + H a + Hire + H a7T . (18) 

The parts contain only electron operators Ci Pr \ for 
the p r orbitals, H a contains only operators for the or- 
bitals Pt,Pzi s i and t ne parts H niJ (H aw ) contain only 

the combinations c] „ rX c%n\ { c \n\ c Wr\ ) witri H = Pt,Pz,s 
running only over the a orbitals. We will use the sym- 
bols Hn, H a , Hira, H an also for the representation of H 
in first quantization, i.e., for complex matrices whenever 
this notation is more convenient. 

In the following, we discuss the 7r-band Hamiltonian 
and the er-band Hamiltonian H a separately before 
we take into account the na hybridization. For this, we 
transform to a fc-space representation of the electron op- 
erators 

C ^ = ^E eiK ' R * C <MA (19) 

where N is the number of unit cells. The preliminary 
definition of the momentum k = kz + K£t has a momen- 
tum component n along the tube and one component Kt 
around the tube. Remember that the alignment of the 
corresponding unit vectors z,t relative to the lattice de- 
pends on the chirality of the tube (see Fig. [T]). 



A. 7r band hopping Hamiltonian 

We start the discussion of the ir band with the limit 
R — > oo, in which the hopping amplitude between all 
nearest-neighbor p r orbitals is V pp . The transformation 
of the hopping terms of the tt band Hamiltonian H£ op to 
fc-space gives 

#hop = V; p J2 w ( K ) c L Pr x^B Pr x + H.c, (20) 

where w(k) = e lK ' Ll + e IK ' L2 + e IK ' L3 . The spectrum 
±V^|io(/e)| of -ff^op i s zer0 at tne two Dirac points 
K = 47r(tcos# + zsin#)/3a and K' = — K. Since we are 
interested mainly in the low-energy states, we expand n 
about each of the two Dirac points, 

K = K + k or K = K' + k, (21) 

to linear order in k = fez + ktt. The resulting approxi- 
mated hopping Hamiltonian for the ir band reads 

H ho P - hv F elT9 ( Tk t ~ ik )cL A T Pr \ c ^BTp T \ + H.c, 

(22) 

with hvp = v / 3|V p p|/2a. The index r = ±1 labels the two 
valleys K and K', respectively. Finally, to bring ([22]) to 
a more convenient form, we change the phase of all p r 
orbitals on the A sublattice by 

CkAr Pr \ -> Te lT0 c kATprX (23) 

and arrive at the usual first quantized form of the Dirac 
Hamiltonian 

H£ ov = hv F {k t a l +kT<j 2 ), (24) 

where the Pauli matrices ui operate in the A,B sublattice 
space. (73 equals 1 on the A sublattice and —1 on the B 
sublattice. Note that, due to the finite circumference 
of carbon nanotubes, k t = (n — tS/3)/R is quantized. 
6 = (iVi - N 2 ) mod 3, for a (N u N 2 )-CNT, is zero for 
metallic nanotubes, to which we restrict the discussion 
in this paper. 

The deviations of t^ Pr from VZ, due to the finite radius 
R leads to an additional contribution to the non-diagonal 
part of Hit, which is added directly to t = V pp w(k), 

{^Wp + V P V e ~ 3lT9 + e3lT9 ( V PP ^ V pv))- ( 25 ) 

This results in a shift of the wave vectors dependent on 
the chirality and radius along the both axes. The curva- 
ture effects also lead to a rcnormalization of the Fermi 
velocity, which can be obtained if we keep terms up to 
second order in a/R and first order in k. The correction 
to fff op is 

AH£ op = h(Av^k t cn + rAvlka 2 ), (26) 
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FIG. 5. (Color online) Curvature induced a bonding between 
p r orbitals of neighbouring atoms (The corresponding atoms 
are defined as A, Bi, and B2). Blue lines are aligned with 
the corresponding p r orbital directions and red (green) lines 
are perpendicular to the red (green) hexagon segment and, 
therefore, to the red (green) bond between the atoms A and 
B2(Bi). Note that all these lines (blue, red, and green) are 
aligned in graphene. Thus, in graphene, only n bonding (that 
is due to overlap of orbitals aligned perpendicular to the con- 
necting line between the two atoms) is possible between such 
orbitals. NT curvature breaks such a symmetry: in curved 
graphene sheet, p r orbital direction (shown by blue lines) has 
perpendicular and parallel components to the line connecting 
two atoms. The latter component is responsible to a a bond 
formation (that is due to overlap of orbitals along the line 
connecting the two atoms). 



with 



HAvp 



2(7V P l + 5V£) + cos66(V; p - V pp , 



128VS 



\R 



hAvp = 



2(VZ, + 3V£) - cos 60(V* - VSL) / a 



128V3 



R 



(27) 



(28) 



Moreover, there is an additional term which couples 
transverse and longitudinal movements, 



AH?_ Z = HAvp- 



with 



hAvp~ 



128^3 Vi? 



((kai + Tk t a 2 ), 



sm6e(v;;-o 



(29) 



(30) 



It should be emphasized that the corrections to the 
Hamiltonian given by Eqs. (|25 p -(|29 p are not the only 
terms responsible for the shifts of the wave vectors and 
the Fermi velocity renormalization. Further terms lead- 
ing to effects on the same order of magnitude indeed arise 
from the tt-ct band coupling. Such corrections by the 
higher energy bands will be discussed in the next sec- 
tions. Collecting all terms we obtain the correction to 



the Fermi velocity 



v F 



a \ ■ 



(31) 



From Eq. ([28]) we find 7 = —0.05, while the direct nu- 
merical solution of Eq. (flj gives 7 = —0.15, which indeed 
reflects the significance of the correction by the energet- 
ically higher bands. In any case, the curvature-induced 
renormalization of the Fermi velocity is small and does 
not give rise to any significant effect in the regime dis- 
cussed in this work. 



B. a band hopping Hamiltonian 

The orbitals s,pt,p z on the sublattices A and B give 
rise to six a bands. For a flat honeycomb lattice (i.e., 
R 00), there is no hybridization between these a bands 
and the it band. For CNTs with finite radii R, finite 
hopping matrix elements between the a and the n or- 
bitals cause a small hybridization which will be treated 
in perturbation theory below. Indeed, one of the small 
parameters appearing in the perturbation theory is the 
ratio a/R with a the lattice constant. It is important 
to understand that, in order to correctly account for the 
curvature effects to order (a/R) 2 in the analytical theory 
derived in Sec. IIV1 it is sufficient to treat the a bands 
in the limit R — > 00, i.e., in zeroth order in a/R. Any 
curvature effect in the a bands of order of a/R would 
lead to a correction of at least third order. 

Furthermore, since the effective theory will be valid 
only near the Dirac points K, K', wc only consider the 
cigenstates and eigenenergies of the a bands at the Dirac 
points, i.e., for k = K, K'. In k space and in the basis 

S={\s A ), \pf), \pf), \s B ), bf), bf)}, (32) 

the a band Hamiltonian has the block form 



H„ = 



where h a ab is given by 



l a,AB 



,AB 



(33) 




2iV s 



-3ir6 



sp 



2rV f 



(V n 

\ r pp 



■ V a ) 

r pp/ 

V a ) 

pp/ 



-3ir0 



-3ir8 



\ pp PP/ 



(V* 

\ pp 



pp) 
V a ), 

r pp / / 



(34) 

r = ±1 labels the Dirac point K,K', respectively H a 
can be diagonalized analytically. We find the energy 
spectrum of the a bands at K, K' 



M*) = ±±ry* -y°\ 
pp 2 ^ pp/* 



" spl 

M*) 

-spl 



1 

2 

spl ' 



± 



^/( es )2 + 18(Vsp) 2 ), 



(35) 
(36) 
(37) 
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The corresponding eigenvectors are given by 



nian H is of the form 



1 



\a pp )^—(-re-^ e \l A '-) + r' + )), 



rB,+\ 



V2 
1 



w spl ) = -T^\ s A )+v+\i B n, 

Ws P 2)=V-\s B )+TV + \l A > + ), 

W* apl }=r V+ \s A )+ V .\l B '-), 



' sp2l 



rj + \s 



T7]-\l 



A,+\ 



(38) 

(39) 

(40) 
(41) 
(42) 
(43) 



where 



\l j ' + ) = ^(\pi)+ir\p{)), 
\l j '-) = j 5 (H)-iT\pi)), 



x/2 



1 ± 



V(£a) 2 + 18(V sp )2 



(44) 

(45) 
(46) 



C. 7ra hybridization 



In this section we consider different mechanisms which 
lead to a coupling between tt and a bands. All these cou- 
plings are weak compared to the typical energy scales of 
the 7r and a bands, so that we may treat them in per- 
turbation theory. It is therefore convenient to write the 
Hamiltonian of a carbon nanotube in first quantization 
in the block form 



H 



(47) 



where is a 2 x 6 matrix describing the na mixing. It 
contains three types of contributions, namely from the 
spin-orbit interaction, from the curvature of the nan- 
otube and from the external electric field, applied to the 
nanotube. As we want to describe the spin-orbit cou- 
pling, we need to explicitly take into account the spin 
of the electron. Thus, each of the 12 matrix elements of 
H-xu is an operator acting on the electronic spin. In the 
matrices describing curvature and electric field effects, 
this operator will be the identity, but for the spin-orbit 
interaction matrices, it will be composed of the spin Pauli 
matrices . Also, the matrix elements may contain mo- 
mentum and position operators. 

In the basis V x S with V = {\p A ), \p B )} and S = 
{\s A ), \pf), \p A ),\s B ), \pf), \ P f)\ the total Hamilto- 



t 
t* 



yAA yAA yAA yAB yAB yAB 

PrS PrPt PrP- PrS PrPt PrP- 

yBA yBA yBA yBB yBB yBB 

p r S PrPt PrPz PrS PrPt PrPz 

7 S 0~ H AB H AB H AB 

n n — h ab h ab h Ab 

U U Sp t PtPt PtPz 

n TjAB ttAB tjAB 

e s 







Vl_l I u I J 

(48) 

where, as explained above, each matrix element is an op- 
erator acting on the electron spin and, furthermore, each 
matrix element may depend on the crystal momentum k 
or on the spatial position r. It will turn out that the only 
spatial dependencies we need to deal with are dependen- 
cies on the azimuthal angle ifi of the carbon atoms. 

The diagonal blocks correspond to the Hamiltonians of 
the isolated n and a bands , as discussed in the two pre- 
ceding subsections, and the off-diagonal blocks (framed 
in bold boxes) are the H vr , and H aTr = H\ a matrices 
with entries V^ v = (w 4 |i?hybK ,J ) for i 7 j <G {A, B} and 
u,v =S {s,pt,p z ,Pr})- Here i?hyb is the part of the 
Hamiltonian inducing transitions between the tt and a 
bands. As stated above, it is composed of three different 
contributions: the hybridization coming from the curva- 
ture H^™ v , from the spin-orbit H^® , and from an ap- 
plied electric field H^ a . Each of these contributions will 
be discussed in the following. 



1. Curvature induced iv-a bond hybridization in nanotubes 

Due to the curvature of the nanotube surface, hoppings 
between orbitals on neighboring carbon atoms, which are 
symmetry-forbidden in flat graphene, become allowed in 
CNTs. This is illustrated in Fig. \5\ We now calculate 
the corresponding hopping matrix elements in fc-space at 
the Dirac points (k = K, K') between the tt orbital (p r ) 
and the a orbitals (s,p t ,p z ). As the hopping is spin- 
independent, the hopping matrix is the identity in spin 
space. Furthermore, since the hopping does not depend 
on the spatial position of the electron but only on its mo- 
mentum, the matrix elements will generally be functions 
of the momentum k but not of the position operator. 
Since we are interested only in the physics near the Dirac 
point, it will be sufficient to consider the Dirac momenta 
k = K, K' only. The neglect of the momentum deviation 
k = k — K^') from the Dirac momenta in the ira coupling 
matrices turns out to be a good approximation; it only 
gives rise to a small renormalization of the Fermi veloc- 
ity, k is only important in the V subspace, as discussed 
in Sec. IHTA"! 

It is important to note that, for a curved surface, 
the vectors L ra connecting nearest neighbors, not only 
have tangential components proportional to z and t, but 
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also a radial contribution L™f. The magnitude L™ is 
of order a/R and thus vanishes in the graphene limit 
R — > oo. The hopping amplitudes must be calcu- 
lated as described in Eqs. (|5j)- (fTU)) , and using the above 
described L n vectors with radial contributions. For the 
transformation to fc-space, on the other hand, the ra- 
dial contributions of L" are not needed because k and 
k arc defined in the two-dimensional tangent space, de- 
scribing the longitudinal (z) and the circumferential or 
transverse (t) direction. Nevertheless, it is convenient 
to write scalar products between the three-component 
vectors L™ = L"z + U/t + L"f and the momentum 

K = KZ + K t t + Of. 

For a given k, the curvature-induced ircr coupling is 
defined by 



rrcurv \ 1 

TV (7 = 2^ 



n=l,2,3 



-i K .L".B„A t 



(49) 



where /i — s,p t ,p z runs over the orbitals forming the a 
band. As we are finally interested in a low-energy theory 
for the 7r band, taking the coupling to the a band into 
accound in up to second order in the small parameters, 
one of which is a/R, and since H^ v will only enter in 
second order perturbation theory, it is sufficient to keep 
only the linear a/R order in H^™ v . Doing so, we find in 
the Px5 basis 



m 



\/3a 
16R 



( 



-2re 3iTe 



V sp 

5VZ,) 

ppj 



_2 Te -3iT9 Vsp \ 

i(3V° p + bv; p ) 

\ r pp r pp I 








(50) 



Again, r = ±1 labels the Dirac points k = K, K', re- 
spectively. Furthermore, note that H™ lv is the identity 
in spin space, so that the matrix in Eq. (|50[) enters the 
block Hamiltonian twice if the electron spin is taken into 
account. 



2. Spin-orbit coupling in ty bands 

For the analysis of the spin-orbit interaction it is im- 
portant to note that the corresponding Hamiltonian H^o 
[Eq. (jf 3[1 ] is local, i.e., it has no matrix elements connect- 
ing orbitals from different lattice sites. Thus, H^® can 
be represented as a 6 x 2 x 2 tensor, corresponding to the 
six a orbitals per unit cell, the two 7r orbitals per unit 
cell and the two-dimensional spin space. Each matrix 
element, however, may depend on the unit cell coordi- 
nate. It is most convenient to write H^® in a local spin 
basis with spin-components S r pointing in the radial di- 
rection of the nanotube, 5* pointing in the transverse 



(or circumferential) direction, and S z pointing along the 
tube axis. Since the local environment, i.e., the defini- 
tion of the direction of the p orbitals, are equal for all 
lattice sites of the tube, the na coupling matrix due to 
spin-orbit interaction has a simple form 



so 



iA 



so 



( 


^ 


s z 





-s l 














s z 


V 





(51) 



Note that there are no S r operators in the matrix el- 
ements of the SOI between the p r and the a orbitals. 
This is easily seen from the second quantized form of the 
atomic spin-orbit Hamiltonian Eq. (|13[) . in which the 
Levi Civita symbol forbids the terms in which a spin op- 
erator and an orbital in the same direction appear. As 
a 7T<7 coupling always involves a p r orbital, the S r spin 
cannot appear. 

In Eq. (fSTj) the directions of the spin operators are 
defined in the local basis f,t,z. On the other hand, the 
hopping Hamiltonian does not affect the electron spin. A 
electron with spin pointing in the global x direction, say, 
and is hopping around the nanotube has its spin pointing 
in the global x direction independently of its position. 
This means that the hopping Hamiltonian is the identity 
in spin space, but only if the global spin basis S x ' v ' z is 
used, rather than the local spin basis S r ' t,z . Thus, also 
the spin operators in Eq. (|5ip must be transformed to 
the global basis. This transformation is given by 



S r (ifi ) = S x cos ip t + S v sin ip t 
S t (ipi) = S v cos ifi - S x sin Lp t 

s z = s z . 



(52) 
(53) 
(54) 



Note that, unlike the hopping matrix H™ lv , ff^o 
pends not on the momentum but rather on ifi . This is be- 
cause the spin-orbit Hamiltonian wc started from did not 
involve hoppings between neighboring carbon atoms but 
only local terms. Transforming this spatially dependent 
part of the Hamiltonian to the momentum space leads to 
non-diagonal matrix elements in k, coupling the trans- 
verse momentum k t to its neighboring momenta kt ± 
As we are finally interested in a low-energy theory for 
the lowest it subband, all virtual +Afc t processes must 
be compensated by a — Ak t process in second order per- 
turbation theory. This will be discussed in detail in Sec. 
TVl For now, wc keep the real-space notation of H^® and 
emphasize again that the matrices in Eqs. (|50|) and (|5Tj) 
are defined with respect to a different basis (fc space and 
real space, respectively). 



Electric field 



As discussed in Sec. Ill C[ there are two significant ef- 
fects of an electric field applied perpendicular to the car- 
bon nanotube. The orbitally diagonal cosine potential, 
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described by H E \ will be discussed in Sec. [V] it turns 

out that H E is reduced by screening effects and, apart 
from leading to a small renormalization of the Fermi 
velocity, has no significant effect on the low-energy the- 
ory we aim at. 

Most important, however, is the s-p r transition, de- 

(2) 

scribed by H E . In the basis V x S, we find for the 
contribution of the s-p r transition to the ira coupling 







\ 







-eE r ((pi)£ 





(55) 



where E r (ipi) is the radial component of the electric field. 
The angular dependence of the electric field in the linear 
response regime can be approximated by 



E r ((fii) = E cos ifii, 



(56) 



where E is the magnitude of the applied electric field 
applied perpendicular to the CNT axis (see also Fig. 1 
in Ref. I60I ). Again, we note that the matrix elements 
of H^ are identities in spin space and do not depend 
on the momentum but on the azimuthal angle ifi of the 
carbon atoms. 



IV. EFFECTIVE HAMILTONIAN FOR THE vr 
BAND 

Having defined all parts of the microscopic Hamilto- 
nian, we are now in a position to derive the effective low- 
energy Hamiltonian of the CNT. This will be done in 
second order perturbation theory. The small parameters 
in which we expand are: the surface curvature a/R, the 
spin-orbit interaction coupling strength Aso = 6 meV, 
and the electric field strength eE£ — 50 • £[V/nm] meV. 
The last two quantities are energy scales and must be 
compared to the typical a band eigen energies which are 
on the order of a few eV. In the derivation of the low- 
energy theory, we completely neglect the Hamiltonian 
H E . For the it band this Hamiltonian alone is known 
to give rise to a Fermi- velocity renormalization of second 
order in the electric field.— The question whether H E ^ 
can actually be neglected will be critically discussed in 
Sec. EH 

The effective Hamiltonian for the it band is calculated 
in second-order perturbation theory as 



I 



Ha 



-Hi 



where 



Tjcurv 



H 



SO 



Hi 



(57) 



(58) 



We proceed by inserting into Eq. (|57[) the unitary matrix 
U ai which diagonalizes H a and is constructed from the 
eigenvectors given by Eqs. (|38 ]) - (|43|) 



m 



Hk + h^uM- 



i 



H, 



■uMhL 



(59) 



so that the inverse operator {e — H a )~ x reduces to a diag- 
onal sum of the a band eigenvalues. Furthermore, since 
we are interested in energies close to the Dirac point, we 
set e = 0. This reduces the complicated 8x8 Hamilto- 
nian matrix, describing it and a orbitals, to an effective 
2x2 Hamiltonian matrix for the it orbitals only. The 
two-dimensional vector space, the matrix H° s is defined 
in, corresponds to the two sublattices A and B. Note, 
however, that the matrix elements of H% s still contain 
the momentum operator k = (fc t , k), the position opera- 
tor ifi and spin operators. In particular, we find 

Hf(0 h k) = hv F [(h + Afc* v )cri + r(fc + Ak z cv )a 2 ] 
+ a[S z a 1 - TS t (ip l )a 2 ]+Tf3 1 [S z cos 36* - S"(^) sin 30] 
+ TeE r (<p i )&[S t (tp i )a 2 -TS z a 1 ] + 

+ 11 TS r (<p i )a 3 + eE r (ip i )Z 1 , (60) 

with the coefficients 



Afci 



v pp y r pp i v pp j 

^ V v pp y pp } w 1 



- cos 39 
sin W 



7i 



Vze s {v; p + v° p )& S o 

18V s 2 p 

v^V^pAso a 

^(^pp ^pp) ^ 
If A 2 
9V 2 ' 

u v sp 

(VZ + V*) a 



R 



, 2 A so 
6 - 



(61) 
(62) 

(63) 

(64) 

(65) 

(66) 



As we aim at an effective theory for the lowest sub- 
band, we project H^. s (ipi,k) onto the subspace spanned 
by the wave functions of this subband. These wave func- 
tions are plane waves oc exp(ikz) which do not depend 
on (p. Thus, projection means (^-averaging the Hamilto- 
nian and setting k t = 0. Furthermore, since there is no 
operator in 7J° ff ((^i,k) which does not commute with fc, 
we may write 

Hf(k) = J ^Hffa(0,k)). (67) 

Any term in -ff° ff (ip, (0, k)) containing odd powers of sin 
or cos average to zero. Finally, the resulting effective 
low-energy theory is given by 

H t = H Lp +#orb + H SO +HsO- ( 68 ) 
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i?orb describes the curvature induced fc-shift of the 
Dirac pointSj 24 i 27 ' 29 For non-armchair metallic CNTs, 
AfcJ v ^ leads to the opening of a gap at the Dirac 
point 



with 



HvpAki 



5.4 meV 

-T^r, ttt COS 30. 



i?[nm] 2 



(69) 



(70) 



The fc-shift Afc^ v along the CNT is irrelevant and has 
been neglected here (see also Ref. IrioT ). 

-ffgQ contains the curvature induced SO I 27 ' 28 which 
does not average out in the ip integration, 



H§& = aS z ax + t(3S z 



with the parameters 



-0.08 meV 
i?[nm] 

-0.31 meV 
Mnml 



cos '39. 



(71) 

(72) 
(73) 



This term contains only on the S z spin operator which is 
consistent with the rotational symmetry of the nanotubc. 

is responsible for the breaking of the electron-hole 
symmetry^ 

Finally, the electric field induced SOI depends on the 
product E{ip)S t {ip), which does not average out in the ip 
integral. Indeed, E(ip) cx cos^ while 5'(</?) = S y cosf — 
S x sin if. The cos 2 if integral leads to a nonvanishing 



H$ = TeEZS v a 2 , 



where 



£ = 6/2 



^so 



3V f 



£o - 2 x 10" 



(74) 



(75) 



Sp 



This term breaks rotational invariance of the CNT and 
involves the S v operator for the electric field along the x 
axis. The fact that the directions of the applied field and 
the induced spin polarization arc perpendicular is typical 
for the electric field-induced SOI. 



V. FIELD SCREENING 

An electric field applied perpendicularly to a carbon 
nanotubc induces a rearrangement of the electrons on the 
tube surface. If the electrons were free to move, as in the 
case of a metal cylinder, the field inside the tube would be 
perfectly screened. However, since the density of states 
in CNTs at half-filling is small, the electrons rather be- 
have as in half-metals, and thus the field screening is only 
partial. In the following, we calculate the field screening 
and its consequences explicitly by two methods: by a lin- 
ear response calculation and by a direct diagonalization 
of the full Hamiltonian, including the external field. The 



analytical linear response calculation will provide us with 
a simple and intuitive picture of the screening but some 
uncontrolled approximations such as, for instance, the re- 
striction to the 7r band, are needed there. The numerical 
calculation takes into account all bands derived from the 
second shell orbitals of the carbon atoms and provides 
a quantitative result, corroborating the linear response 
calculation. 

The task is, therefore, to calculate the charge response 
of a CNT to the electrostatic potential described by 



He = 0o,tot cos(<y9 niC )n niC 



(76) 



with </>o,tot the amplitude of the total electrostatic po- 
tential acting on the tube surface. <p n ,C IS the azimuthal 
angle of the carbon atom in unit cell n and sublattice Q. 
The operator n n x = 2u c n c M Cn ^ ! ^ counts the electrons 
on atom (n, £) in all second shell orbitals fi = s,p r ,p t ,p z . 
In order to keep the notation simple, we drop the spin 
index A in this section and multiply the charge response 
by a factor of 2. 

Note that He describes a homogeneous electric field E, 
which gives rise to the electrostatic potential cx t(<£>) = 
eER cos if. E is the field which is applied externally. 
In the following we will show that the induced electron 
charges on the tube surface give rise to an induced elec- 
trostatic potential 0i n( j (r) which, although it has a rather 
complicated spatial structure away from the nanotubc 
surface (i.e. for |r| ^ R), reduces to ^ind(v) = <A),ind cos <f 
for |r| = R. Thus, anticipating this result, the total 
potential 0tot(¥>) = ^ext(y) + 0ind(<^) = ^o.tot cos is 
described by the Hamiltonian (f76")) . As discussed in Sec. 
IV Dl this additivity property of the amplitudes 0o,ext and 
</>o,ind can be interpreted as a linear screening of the elec- 
tric field inside the nanotube, i.e., 0o,tot = eE*R with 
the screened field E* < E. 



A. Linear response of the tt band 

We define the static linear response coefficient 
Xqqi (n, n') as the proportionality constant between the 
density response at site (n, £) in orbital /i and the density 
perturbation at site (n', £') in orbital // described by any 
Hamiltonian -Ac^,,, „,c n '^ ^ 

" t (77) 



P„0 = 6 [cl^Cn^j = Ax££ (n, n'). 
X(£i (n, n') can be calculated by the Kubo formula 

X $'(n J n , ) = -2i 



die - * 



(78) 



JtH -itH 

rcscntation of the electron annihilation operator and 



where c n ^ iM (t) = e ltH c n ^ jAI e %tti is the Hciscnbcrg rep- 
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r/ = + ensures the convergence of the time integral. The 
factor 2 accounts for the spin-degeneracy. The average 
denotes the expectation value with respect to the ground 
state of the electronic system. The generalization to fi- 
nite temperatures is possible but not of interest here. 

The linear response of the nanotube to the Hamilto- 
nian defined in Eq. ()76[) is 



E 



cos(<p n ',c')Xcc' (n,n') 



(79) 



In the remainder of this subsection we deal only with 
the linear response in the n band allowing us to set 
(! = (/ = p r and suppress the orbital index. The con- 
tributions of the a band will be discussed below on the 
basis of the numerical calculation. Furthermore, we re- 
strict the calculation to armchair nanotubes. 

A straightforward calculation of the charge density in- 
duced by the Hamiltonian (f76"]) gives rise to two terms in 
the charge response: a normal response, which has the 
same cosine modulation as the inducing Hamiltonian [Eq. 
([75)) ] , and an anomalous response, which is staggered on 
the sublattice level and has a sine modulation, 

P(<Pn,t) = <fo,tot [Xn COS - (Xa,Bm<p n ,c] , (80) 

with the normal and anomalous response coefficients 

Xn = Xaa + Rcxsa cos(ip AB ) - Im^sA sin(ip AB ), (81) 
Xa = Rcxba sin(^As) + ImxBA cos(ipab)- (82) 

Here fAB = af^/SR is the difference of the azimuthal 
angle between the different sublattice sites within one 
unit cell and xaa,Xba are the fc-space susceptibilities, 
defined by 



x£? r (q = t/-R)- 



(83) 



As expected, the relevant susceptibilities for the charge 
response to a transverse homogeneous field are to be eval- 
uated for zero momentum along the tube and for the 
smallest possible momentum around the tube. 

The charge susceptibilities required for Eq. ([53")) can 
be expressed as fc-space integrals 



xcc 



Q(e F - e (/c))8(e /(K + q) - e F ) 



e a (n) - e a >(K + q) 



(84) 



-cc(q) = a «' ex p«(^ - i+ q )) 



x Q(e F ~ e a (K))®(e a ,(K + q) - e F ) 
e a {K) - £ a j(n + q) 

where the spin-degeneracy has been taken into account. 
£p is the Fermi energy and e a ('*) are the energy eigen- 
values for the two branches (a = ±1) of the tt band. 
We assume here that the curvature effects, discussed in 
Sec. IIII1 do not affect the final charge polarization signifi- 
cantly. This assumption will be tested numerically in the 



next subsection. Thus, we calculate e±(re) = ±|V^,u;(k)| 
and 4> K = argw(/c) from the Hamiltonian ([2"U)| . 

Furthermore, we restrict the discussion to the charge 
neutrality point e F = so that a = — 1 and a' = 1 and 



x-cc(q) 



\V^\N ^ K«)| + KK + q) 

1 ^ eX P (^(^K - <?Wq)) 



\YS,\N*-< \w(K)\ + \w(K + q) 



(87) 



Note that for q = we have Xcc(0) = ~X-Cc(0) so that 
the linear response to a homogeneous potential is zero in 
flat graphene. This is a consequence of the vanishing den- 
sity of states of graphene at the charge neutrality point. 
For the q vectors given by Eq. ([53")) and for the fc-space 
grid defined by the circumference and the length of the 
CNT (we assume periodic boundary conditions in z direc- 
tion), the fc-space summations in Eqs. ((86)) and ([87)) are 
evaluated numerically. For instance, for a (10,10)-CNT 
and in the limit of infinitely long tubes we obtain 



\VZ\Xr. 



-0.0424, 



\VIJxa = 0.002. 



In the Dirac approximation of the band structure, i.e., 
e±(k) = ±V3|V^||ka|/2, Eqs. (JSBJ) and (J57J) can be 
calculated analytically for small q — |q|. For the spin- 
susceptibility, which is equal to the charge susceptibility 
for non-interacting systems, this has been done in Rcfs. 
I62T4641 with the result 

Hvp , s Aa qa 

+ 16 
Aa 3qa 

2^ ~ ~16~' 

where A is an ultraviolet cutoff, a is the lattice constant, 
and hvp/a = \/3|V^p|/2. From this we find the analytical 
forms for the normal and anomalous response coefficients 

\V P l\ Xn = -^-^(^Y + 0((a/Rf) (91) 



a 

hvp 



-xcdi) 
x-tdv) 



(89) 
(90) 



iV^R 6\/37r \R< 

ra^ = sl4(s) 2+0 (w fl ) a ) 



(92) 



The comparison of the analytical result for Xn with the 
numerical evaluation of the Kubo integrals in Fig. [5] 
shows that the first term in Eq. ([9"T)l captures the lead- 
ing a/i?-term qualitatively in a correct way. However, 
the Dirac approximation leads to a wrong prefactor. 

The analytical form of the anomalous response calcu- 
lated in the Dirac approximation differs from the exact 
result. In the next subsection we explain the reason for 
this and why the evaluation of Kubo formulas within the 
Dirac approximation is not reliable. However, as we are 
finally interested in the field screening, the anomalous 
response, which is staggered on the atomic scale, is irrel- 
evant as it averages out in a continuum approximation 
of the charge distribution. The staggered potential in- 
duced by the anomalous response has a sin ip modulation 
so that it does not open a gap. We have checked numer- 
ically that such an additional term in the Hamiltonian 
does not change our results. 
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FIG. 6. (Color online) The normal response coefficient Xn of 
the 7r band as a function of a/R. The black dots are the exact 
result from the numerical evaluation of Eqs. (|86|) and (|87|) . 
The solid line shows the linear term in Eq. ()9ip . The inset 
shows the anomalous response coefficient \a together with a 
polynomial fit \a — 0.016(a/i?) 2 . 



B. Analysis of the Kubo integral 

It is instructive to study the structure of the Kubo in- 
tegrals for the charge susceptibilities [Eqs. (jSo) and ([57])] 
in more detail, especially in view of their Dirac approxi- 
mations. 

For q = the intra-sublattice susceptibility has the 
structure 



xcc(q = °) K J de 



D(s) 



(93) 



where D(e) = TV 1 J2 Ka ^( £ ~ £ a{^)) is the density of 
states. In graphene, D(e) is known to be linear in the 
energy e near the Dirac point. For carbon nanotubes, 
which can be viewed as graphene with one confined di- 
rection, D(e) is constant for e ~ 0, but it increases dis- 
continuously with ease crosses more and more transverse 
subbands. On a coarse grained energy scale D(e) ~ e for 
both, CNTs and graphene, as long as e < 3eV. In the 
Dirac approximation, it is assumed that only electronic 
states with low energies are important, i.e., that the inte- 
gral in Eq. (|§3"|) converges before e ~ 3eV. However, Eq. 
(|93p docs obviously not converge for D{e) ~ e so that, 
strictly speaking, the Dirac approximation is not allowed 
for Kubo formulas. 

This convergence problem is reflected in the cutoff de- 
pendence of the first terms in Eqs. (|59")l and How- 
ever, the terms linear in q do not depend on the ultravi- 
olet cutoff, which suggests that only contributions from 
small energies, where the Dirac approximation is valid, 
enter the g-dcpcndcncc. Indeed, one finds that in the 
high energy regime 



xcc(q) - xcc(°) 

converges quickly. 



D( 



(94) 



In order to compare the Kubo integrals of the tight- 
binding formulation with the Dirac approximation, we 
consider xaa(<\ = 0). Eq. (|86|) can be written as 



(95) 



|V£|X«(0) = / de/(e 
Jo 



with 



f(e) 



TV ^ 



S(s-\w(k)\) 
2\w(k)\ ' 



(96) 



In the Dirac approximation in which w(K + k) oc k x + 
iky is assumed (we may drop all the prefactors in this 
analysis in favor of notational simplicity), /Dirac (e) can 
be calculated easily. Since /(e) ~ D(e)/e, the Dirac 
approximation /Dirac (e) must be constant. 
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FIG. 7. (Color online) The energy-resolved integrand of the 
Kubo integral /(e) for Xcc( < l ~ 0)- The black line is the inte- 
grand calculated from the full tight-binding model /tb(e) (see 
text). The horizontal gray (red) line shows the integrand in 
Dirac approximation /Dirac (e). The inset shows the integrand 
g(e) for q = 27rt/50. The dots are the results of the tight- 
binding calculation for a 200 x 200 fc-space grid and the line 
shows the Dirac approximation of g(e). 

Calculating /(e) numerically for the more complicated 
tight-binding form of w(k), which can be evaluated only 
for finite size systems TV < oo, is not quite straightfor- 
ward. We do this by discretizing the energy in Ae steps 
and defining 



ftb(e,N) 
1 



TVAe 



E 



Q[\w{k)\ -e]9[e + Ae - \w(k) 
2K«)| 



(97) 



In the limit Ae — > and TV — > oo, Eq. (|97|) approaches 
the actual /(e), as defined in Eq. However, finite- 

size effects make f t b look very rugged. In order to obtain 
a smooth /(e), which can be plotted nicely, we define an 
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average over N 



7r band only 



1 M 
/tb(e) = E 



/ tb (e, (202 + 3m) 2 ) + / tb (e, (203 + 3m) 2 ) 



(98) 



where we have assumed a quadratic fc-space grid and we 
ommitcd all grids in which a fc-space point hits a Dirac 
point. The choice of the smallest N = 202 2 is completely 
arbitrary, as well as the choice of M = 82, as long as M 
is large enough. 

Fig. [7] shows /Dirac (e) and /tb(e)- Obviously, for small 
e <C \Vp p \ the Dirac approximation coincides with the 
tight-binding calculation. However, for larger energies 
they strongly differ. Moreover, while the finite tight- 
binding bandwidth < 3) gives a natural high en- 
ergy cutoff and ensures the convergence of the energy 
integral, the Dirac model must be replenished by an ul- 
traviolet cutoff A in order to ensure the convergence. 

For the deviation of the finite q susceptibility from 



Kp|[xcc(q)-xcc(q = o)] 



deg(e), (99) 



with 



\w{k)\) 
1 



\w(k)\ + \w(k + q)\ 2\w(k)\ 



, (100) 



the integral is convergent for the tight-binding formula- 
tion as well as for the Dirac approximation. Furthermore, 
17(e) is only large at low energies, i.e., where the Dirac ap- 
proximation is valid. Also, the comparison of the Dirac 
approximation of g(e) and the tight-binding calculation 
in the inset of Fig. [7] shows that the integrand g(e) is 
equal in both calculations. 



C. Exact charge response 

In order to scrutinize the results of the linear response 
calculation we diagonalize the Hamiltonian Hq + He nu- 
merically and calculate the charge distribution induced 
by He [Eq. ([To]) ]. Unlike the linear response method, 
this calculation is not restricted to small fields E. We 
start with considering only the 7r bands in order to be 
able to check the results of the linear response calcula- 
tion directly. In a second step, we then take into account 
all carbon orbitals of the second shell in order to see how 
the charge response is affected by the a orbitals. 



Considering only the ir band and neglecting all cur- 
vature effects, we set H = H£ . We consider only 
armchair nanotubes and transform the direction along 
the nanotube to fc-space. For a nanotube with N c unit 
cells in circumferential direction, the Hamiltonian is a 
fc-dependent 2N C x 2N C matrix, which wc diagonalize nu- 
merically. From the eigenvalues e m> k and the correspond- 
ing eigenvectors ipm^kif) we calculate the induced charge 
density (in units of the electron charge e) 



£F 



■,k)\i l m,k( i P)? 



1. 



(101) 



Note that <p must be considered as a discrete variable 
with 2N C possible values between and 2ir correspond- 
ing to the A and B sublattice sites in the N c unit cells 
in circumferential direction. N z is the number of unit 
cells along the tube (z direction). The spin-degeneracy 
is taken into account in Eq. (|101j) . 




FIG. 8. (Color online) The induced -k band charge density 
p(ip) for in a (10,10)-CNT for o ,tot = 0.01 V^. The black 
dots show the results of the numerical diagonalization of the 
7r band Hamiltonian. The line (red) is the charge response 
evaluated in linear response. The linear response data is also 
discrete. It is joined that it can be distinguished from the 
numerical results. Actually both calculations give the same 
result. 



Fig. [8] compares the density response to an external 
field corresponding to (/>o,tot = O.OlVpp, calculated from 
linear response and from the direct solution of the lat- 
tice Hamiltonian. Within numerical accuracy, both cal- 
culations give the same results for these small potential 
amplitudes. Also, one can clearly see the anomalous re- 
sponse, which makes the response curve different from a 
pure cosine shape. 

The response coefficients Xn and Xa can be extracted 
from the numerical calculations of the induced charge 
densities p(<p) by a fit to Eq. ([50)) . The normal response 
coefficients Xn f° r larger potentials </>o,tot resulting from 
these fits are shown in Fig. [S] For small external po- 
tentials, Xn(0o.tot) is close to its linear response result, 
Eq. (|9"Tj) . For larger potentials, however, there are large 
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0o/A sb 

FIG. 9. (Color online) Nonlinearities in the normal response 
coefficient \ n in the -n band. The different curves correspond 
to different nanotube radii R — 0.67, 1.00, and 1.33 nm. The 
solid lines are the response coefficients fitted to the numer- 
ical charge density calculated by Eq. (|101|l . The dashed 
lines indicate the corresponding linear response coefficients. 
The applied potential (abscissa) is rescaled with the subband 
splittings of the nanotubes with the different radii. These are 
A sb = (0.31, 0.16, 0.078) | Vpp | for the radii R = 0.67, 1.00, and 
1.33 nm, respectively. 



deviations. The linear response regime is left if 0o,tot is 
on the order of the subband splitting A s b- This subband 
splitting depends on the CNT radius and is given in the 
caption of Fig. |U 



2. a band contributions 

In the analysis of the Kubo integral we have observed 
that the contributions to the charge susceptibilities from 
high energies cannot a priori be neglected. Thus, in order 
to check the stability of the ir band calculation of the in- 
duced charge density, we calculate the charge redistribu- 
tion in a tight-binding model containing all second shell 
orbitals of the carbon atoms, i.e., the n and the a bands. 
For this, we numerically diagonalizc the Hamiltonian 

H = H bop + H so +H E , (102) 

with the individual terms of H defined in Sec. |TTJ The 
charge density can be separated into a n band part and 
a a band part 

^(^n, C )=E( C <P^nC P ,A) (103) 
A 

/^n, C )=E E ( C nCMA C nOA) • (104) 

Both terms have the qualitative form found in the linear 
response theory. Defining the normal and the anomalous 
response coefficients for the 7r and a response separately, 



we find for the exemplary (10,10)-CNT 

K P \xl = -0-050 V; v xl = -0.0077 (105) 
V; p Xa= 0-011 V; pX : = 0.0049 (106) 

Thus, about 87% of the normal charge response, which 
will be important for the field screening discussed in the 
next subsection, comes from the it electrons. Note that 
the field screening, which enters the Hamiltonian of this 
numerical calculation, has been taken into account self- 
consistcntly. 

D. Field screening 

As mentioned before, the charge response calculated in 
the previous subsections is not the actual charge response 
of the carbon nanotube because electron-electron inter- 
actions have not been taken into account. Their effect, 
however, is important because they reduce the amplitude 
of the charge rearrangement considerably 

In order to see this, we include the electron-electron in- 
teractions with a self-consistent Hartree calculation. In 
addition to this mean-field approximation, we assume 
that the charge induced by an external electric field is 
distributed homogeneously on the carbon nanotube sur- 
face, i.e., we neglect the anomalous response as it is stag- 
gered on the sublattice level and averages out in the con- 
tinuum approximation. We therefore assume that the 
electric field induces the continuous charge distribution 

p 3D (r) = p cos(ip)S(\r\-R), (107) 

where R is the radius of the nanotube, <p is the continu- 
ous azimuthal angle of the spatial coordinate r measured 
from the center of the nanotube, and po is the ampli- 
tude of the induced charge distribution, po is most easily 
determined by requiring that the total induced charge 
in the positive half space (cos ip > 0) calculated from 
the continuum approximation Eq. (|107[) and from the 
lattice-resolved expression Eq. ([50)) are equal, 

/ d 3 r6(cos</?)p 3D (r) = e ^ 6(cos <p n £)p((p n ,c), 

J n,< 

(108) 

with Q(x) the unit step function and e the electron 
charge. The anomalous response drops out on the right 
hand side of Eq. (|108[) . Furthermore, for the normal 
response, the right hand side of Eq. (|108[) is calculated 
approximately by using 

n,C a V J 

where Aip = is the mean azimuthal angle between 
neighboring lattice sites in the armchair CNT, L is the 
length of the tube (in z direction) and f(<p) is any func- 
tion. The consequences of the errors, introduced by the 
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above approximation, for the screening properties are 
much smaller than the consequences of the uncertainties 
of the hopping parameters. Therefore, we have 



Pa 



4e 



■Xn4>0,tot, 



(110) 



with </>o,tot the amplitude of the electrostatic potential 
at the tube surface, including the contributions from the 
external and the induced fields. This back action effect 
is the basis for charge screening and will be calculated in 
the following. 

A charge distribution of the form (|107[) induces an elec- 
tric potential (f>i n d which can be calculated by 



<And(r) 



4tt£ 



(111) 



where £o is the vacuum dielectric constant. For symmetry 
reasons, </>; n d does not depend on the z coordinate (along 
the tube) but has only a radial (r) and azimuthal (tp) 
dependence. We find 

<t>ind(r, <p) = </>o,ind cos(ip)f{r/R). (112) 

with O) ind = -I^o.tot and 

/(^) = {| IVrtn- < 113 > 

The proportionality constant T is given by 

r = 356.78i?[nm]|x„|eV (114) 

The total electrostatic potential 0tot(r) felt by the elec- 
trons on the tube surface consists of two parts 



</>tot(r) = ^ext(r) + ^ind(r) 



(115) 



where ^> e xt(r) comes from the external electric field and 
</'ind(r) is induced by the rearrangement of electron 
charges at the tube surface. 

At the tube surface |r| = R, all terms in Eq. (|115[) 
have the same functional form (~ cos (p), so that the 
problem of solving the self-consistency equation reduces 
to an equation for the amplitudes of the electric poten- 
tials at the tube surface 

00, tot = 0O,oxt + 0O,ind = 0O,cxt — I>o,tot (H6) 

and so 

0O,tot = 0O,oxt/7; 7 = 1 + T. (117) 

Inside the tube (r < R), the functional form of </>tot (r) is 
the same as </f> ext (r), so that the total electric field inside 
the tube is homogeneous, just as the external field, but 
is screened by the factor 7. Thus, we may define the 
screened field 



E* 



E/j. 



(118) 



Note, however, that only inside the nanotube, the field 
is screened homogeneously. Outside the tube, the total 
electric field is highly inhomogeneous. 

The screening factor 7 quantifies how free the surface 
charges arc to move. For a metal cylinder 7 — > 00, which 
means that the charges can move freely on the surface. 
For a cylinder made from an insulating material, the 
charges are localized and can only form dipoles, which 
is reflected by 7 > 1. A carbon nanotube lies in between 
the metallic and the insulating limit. For a (10,10)-CNT, 
for instance, we find 7 ~ 5.6, where the tt and a bands 
are taken into account. Note that 7 — 1 depends lin- 
early on the inverse hopping parameters [see, e.g., Eqs. 
(|105[) and (|106p ]. so that the typical error of 7 in a tight- 
binding calculation as performed here can be estimated 
by the spread of tight-binding parameters found in the 
literature. Within these error bars of 30%, our result is 
in agreement with previous worksi 65 ! 66 



VI. ANALYSIS OF THE LOW-ENERGY 
THEORY 

A. Spectrum without electric field 

In the absence of electric fields, the CNT is invariant 
under rotations around its axis and the z projection of the 
spin s z = ±1 is a good quantum number. In addition, the 
valley index r = ±1 is always a good quantum number, 
as long as there is no intervalley scattering. The term 
hv F Ak z v a2 in H™ b shifts the momentum of the Dirac 
point along the tube. Assuming infinitely long tubes, 
this shift is irrelevant and can be dropped. Thus, we are 
left with the Hamiltonian 

H = Ts z p + Thv F k<T 2 + {hv F (k t + Ak t cv ) + as z )a 1 , (119) 



which is readily diagonalized (see also Ref. 123). The 
eigenvalue spectrum has eight branches, given by 



£u. s - (fc) = ts z (3 + y/{hv F k) 2 + (hvpik* + Afc* v ) + as*) 2 , 

(120) 



e d s (fc) = TS z f3 - y/(hv F k) 2 + {hv F {k l + Afc* v ) + as 2 ) 2 , 

(121) 

where r = ±1 and s z = ±1. In the basis 
{\A t) , \B t) , \A I) , \B I)}, the eigenvectors are given by 



1 



-1 











d+ \ _L —I 

/2 I I ' V2 1 ' V2 1 ' V2 -1 



(122) 



with 



e i# ± = (hvFjk 1 + Afc« v ) ± a) + irhvp k 

^(hvpikt + Afc* v ) ± a) 2 + (hv F k) 2 ' 
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In the case of a semiconductor CNT a gap between the 
two nearest spin-up and spin-down states is 



0.4 



A = 2\a±0\, 



(124) 



where ± corresponds to electrons and holes, respectively. 
In the case of an armchair nanotube, the eigenstates are 
two-fold spin-degenerate and the gap between holes and 
electrons is 2\a\. 

An electric field perpendicular to the tube axis gives 
rise to the spin-orbit term i?so °c S y , and so s z no longer 
is a good quantum number. In the following we discuss 
the consequences of such an electric field. 



B. Helical modes 

As already discussed in Ref. [6(1 the interplay of strong 
electric fields and spin-orbit interaction leads to helical 
modes. For an armchair CNT the chiral angle is 8 = tt/6 
and all terms proportional to cos 36* vanish, i.e., Afc* v = 
and /3 = 0. Furthermore the longitudinal fc-space shift 
Ak z v can be ignored since it can be removed by regaug- 
ing the phase of the orbitals. This leads to the effective 
Hamiltonian for the lowest subband in an armchair CNT 



H &xm = Thv F ka 2 + aS z a x + TeE£S y a 2 , (125) 
which has four branches of eigenvalues 



e(fc) = ±eE£ ± yfa 2 



(hv F ky 



(126) 



for each valley. In the basis {\A f) , \B f) , \A I) , \B I)}, 
the corresponding eigenvectors are given by 




1 -e" 




(127) 



with 



a + irhvFk 
\J o? + (rhvFky 



(128) 



Note that the eigenvectors are independent of the electric 
field E; only the energetical order depends on E. In the 
following, we label the four branches for each valley by 
n=l, 4. For each k, n = 1 corresponds to the highest 
eigenvalue and n — 4 to the lowest. 

In Figs. [TUTa 1 ) and lTTT a). the spectrum calculated from 
the analytical low energy theory is compared to the spec- 
trum calculated from the numerical solution of the com- 
plete tight-binding model [Eq. (JT])]. Obviously, there are 
significant differences. These are due to the neglect of 
the Hamiltonian Hg in the derivation of the analyti- 
cal theory. It is known^i that alone gives rise to a 
renormalization of the Fermi velocity on the order of the 
squared screened electric field (E*) 2 = (E/j) 2 . However, 
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FIG. 10. (Color online) (a) Low-energy spectrum for a 
(10,10)-CNT (armchair) in a field E = 1 V/nm. (b) The 
dependence of j/-spin polarization {3,k\S v |3, fc) on the mo- 
mentum k along the tube. The solid lines are the results of 
the analytical low-energy effective theory. The dots corre- 
spond to numerical calculations (see text). 



the combination of with other parts in the spin-orbit 
Hamiltonian which depend trigonomctrically on ip may 
give rise to effects of first order in E* ~ l/j. Thus, 
the differences of the numerical and analytical spectrum 
should disappear in the limit of perfect screening, i.e. 
7 — > oo. In order to check this expectation, we increase 
the screening parameter 7 in the numerical calculation 
and find that the agreement between numerics and an- 
alytics becomes better for larger screening. For 7 = 20, 
the analytical theory is as good as the numerics, as shown 
in Fig. [T2j Better screening can be achieved, for instance, 
by filling the CNT with a dielectric. 

The reason why these effects arc important although 
they arc of higher than second order is that the parame- 
ter regime we are discussing here is at the border of the 
applicability of perturbation theory. Nevertheless, as we 
will argue now, the effective theory is valuable also in 
this regime because it correctly captures the most im- 
portant features of the helical modes. Moreover, we will 
show that the numerical results predict an even better 
spin-polarization than the analytical theory. 

As mentioned above, the actual value of the parame- 
ter £ is unknown. We use here a conservative estimate 
£ = 2 • 10~ 5 , but exhausting the range of values for £ > 
Aso> etc., £ can increase by one order of magnitude.— In 
this case, i.e. for £ = 2 • 10~ 4 , the analytical theory fits 
much better to the numerical calculations. This is essen- 
tially because the spin-splitting generated by eE£S z <T2 is 
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FIG. 11. (Color online) a) Low-energy spectrum for armchair 
CNT (20,20) in a field E = 1 V/nm. b) The dependence 
of y-spin polarization (3, k\ S y |3, k) on a wave vector k. The 
solid lines are the results of the analytical low-energy effective 
theory. The dots correspond to numerical calculations (see 
text). 




FIG. 12. (Color online) Low-energy spectrum for armchair 
CNT (10,10) in a field E = 1 V/nm assuming that the screen- 
ing factor 7 is equal to 20. The solid lines are the results of 
the analytical low-energy effective theory. The dots corre- 
spond to numerical calculations. In this limit we observe a 
good agreement between two models, which confirms our hy- 
pothesis that discrepancies in Fig. [Tolare caused by the fact 
that we neglected the change in charge distribution caused by 
electric field. 



sufficiently large so that the higher order terms, causing 
the deviations from the simple analytical model, are not 
effective in this case. We provide a movie in the EPAPS 
which illustrates this.— In this movie, the numerical spec- 
trum is compared to the analytical spectrum as the pa- 
rameter £o is varied between 0.5A and 0.05A. It is seen 



that the agreement is better for larger £o- 

Next, we discuss the spin-polarization. The eigen- 
states given in Eq. (|127p allow us to calculate the spin- 
polarization of each branch in the analytical model. It is 
easily seen that (k, n\ S x,z \k, n) = 0, i.e., the spin is only 
polarized in the y direction (perpendicular to the CNT 
axis and to the electric field). In this sense, the spin is 
perfectly polarized, even though (S y ) is smaller than one; 
the vector (S), with S = (S x , S y , S z ), is perfectly parallel 
to the y direction and has no components perpendicular 
to y. This j/-spin polarization in the analytical model is 
given by 



(k,n\S y \k,n) = ± 



VA: 2 + (a/hv F y 



(129) 



Note that (k, n\ S y \k, n) is odd in k. This means that 
for Fermi levels as indicated in Figs. [T0]and[TTl we have 
one helical liquid per Dirac point. The sign of the helic- 
ity, however, is the same for each Dirac point (this is a 
consequence of time-reversal invariance) , so that a CNT 
with a properly tuned Fermi level is a perfect spin filter. 

In Figs. [TOlb) and [TTTb). the analytical spin- 
polarizations are compared with the numerical spin- 
polarizations. They agree well and the analytical spin- 
polarization is seen to be a lower bound to the numerical 
result. One should also note that, because of the in- 
creased splitting at zero k between states 3 and 4 of the 
numerical calculation compared to the analytical results, 
the range of possible Fermi levels for the helical liquid 
is increased. This leads to a higher maximum (S y ) in 
the more rigorous numerical solution. Also in this sense, 
the analytical model provides a lower bound on the max- 
imum (S v ). For example, in the case of a (10,10)-CNT 
to which an electric field E = lV/nm is applied (see Fig. 
fl7j|) one can achieve (S y ) ~ 90% in the helical phase. 

Another well studied effect of the Hamiltonian 
is a rcnormalization of the Fermi velocity We have not 
taken this rcnormalization into account in the analytical 
calculation in Fig. [TTJ but it is accounted for in the 
numerical calculation. Also this effect allows us to go to 
larger k F in the helical regime. 

Now that we have understood which features are 
well captured by the analytical model (e.g., the spin- 
polarization and the spectrum of the holes), and which 
are not described correctly (the electron spectrum for 
large fields), we discuss the case of non- armchair nan- 
otubes on the basis of the analytical model. As men- 
tioned above, the effects leading to deviations from the 
analytical results can be suppressed by increasing the 
screening of the externally applied electric field. 

For non-armchair, but metallic nanotubes (i.e., k t = 
0), the chiral angle is 8 ^ tt/6 and cos(3#) ^ 0. This 
gives rise to two additional terms in the Hamiltonian [see 
Eq. (|125|) ]. One term, Hv F Akl v ai, opens an orbital gap. 
The other term, t/3S z , acts as an effective Zeeman field 
along the tube axis. For the mechanism of valley suppres- 
sion, discussed in the next subsection, it is important to 
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note that both additional terms have opposite signs in 
different valleys K, K'. 

The most prominent effect of 9 7^ tt/6 is the opening of 
a large orbital gap in the meV range for reasonable CNT 
radii R. For a nanotube with R = 1.44 nm but different 
helical angles, the low-energy spectrum is shown in Fig. 
1131 In the armchair CNT, the helical liquid appears on an 
energy scale of a few hundred ficV. Thus, it is not a priori 
clear that the helicity of the left/right movers survives 
the departure from the pure armchair topology of the 
lattice. However, a plausibility argument for the stability 
of the helicity with respect to hv f Ak^ai can be given: 
for large k 3> Afc* v , the electronic states are eigenstates 
of the operator <J2S y . In this limit the y-spin polarization 
(i.e., the helicity) becomes 100%. The y-spin polarization 
changes sign under k — > —k. Furthermore, as long as 
there is no band crossing, (S y ) must be a smooth, odd 
function of k. This means, as long as \kp\ 7^ 0, the 
y-spin polarizations of the second band at ±kp must be 
opposite. The question is only, how large is the amplitude 
of the polarization. 




FIG. 13. (Color online) Low-energy spectrum for two chiral 
metallic CNTs with different chiralities but comparable radii. 
The field strength is 0.5 V/nm. The four bands with smaller 
absolute energy (red, solid) correspond to a (23,20)-CNT with 
R — 1.44nm and 9 = 0.154-7T. The four bands with larger 
absolute energy (blue, dashed) are from a (26,17)-CNT with 
R = 1.44nm and 9 = 0.128tt. 



The orbital gaps are not the only effect of a nontrivial 
chirality 6 7^ ir/6. The effective Zeeman field rj3S z leads 
to an additional spin-polarization in z direction (along 
the tube). This can be observed in Fig. Q3J At small k, 
where the electronic states are not forced into eigenstates 
of (T2, the effect of this effective Zeeman field is largest, 
i.e., the z polarization is maximal, while, at large \k\, 
the spin tends to be aligned in y direction. This effect 
reduces the quality of the spin helicity in that the spin 
alignment is not completely odd in k and not along y. 
Only the y-spin component is odd, but the z-spin is an 
even function of k. 
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FIG. 14. (Color online) Spin-polarizations of the third band 
in y and z direction in chiral nanotubes as functions of k. 
The field strength is 0.5 V/nm. The odd functions around 
k — (red) correspond to {S y ) and the even functions (blue) 
to {S z }. The solid lines are results from a (23,20)-CNT and 
the dashed lines from a (26,17)-CNT. 



C. Valley suppression 

In the previous section we have seen how the quality of 
the helicity is reduced in case of a non-armchair chirality 
of the CNT. This was due to an orbital fc-shift hvp/S.k^ai 
and an effective Zeeman field t/3S z , both of which are 
consequences of the chiral angle 9 ^ tt/6 deviating from 
the armchair case. Here, we show that the appearance of 
these terms can be turned into an advantage: by applying 
an additional magnetic field along the tube, the perfect 
helicity can be restored in one valley while all bands of 
the other valley are removed from the low-energy regime. 

First note that both, the orbital fc-shift and the effec- 
tive Zeeman field, can be interpreted as originating from 
a magnetic field applied along the axis of the CNT. In 
general, a magnetic field B has two effects on the elec- 
trons. First, it induces a Zeeman energy, described by 
the Hamiltonian 



H z = fi B B • S, 



(130) 



where \xb is the Bohr magneton and S is the vector of 
Pauli matrices for the electron spin (eigenvalues ±1). 
Second, if the magnetic field has a component along the 
CNT axis, i.e., B z 7^ 0, the transverse wave function of 
the electron encloses magnetic flux and this gives rise 
to a shift in the electron momentum in circumferential 
direction 



Ak% 



ttB z R 



(131) 



where $0 



Ak* B must be added to Afc' v , 



h/\e\ is the magnetic flux quantum. Since 
as well as fiB z S z must be 
added to t/3S z , the effect of a non-armchair chirality of 
the CNT can be compensated for by a magnetic field, at 
least to a certain degree as explained below. However, 
since the chirality-induccd Afc* v and t/3 have opposite 
signs in different valleys, but the real magnetic field terms 
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Afc^ and /j,bB z have not, this compensation works only 
in one of the two valleys. In the other valley, instead, the 
effect of the chirality is even increased. 

Note that in general only one of the two chirality ef- 
fects can be compensated for by a magnetic field B z . 
Compensating for the orbital gap requires 



ftp. 



nB z R 



5.4meV 
t— — cos SO, 



(132) 



$ ' i?[nm] 2 
while compensating for the effective Zeeman field requires 



0.31mcV 
t — — — cos SO. 



R[i 



(133) 



In general, Eqs. (|132|1 and (|133[) are not compatible. 

However, these two conditions have a different depen- 
dence on the CNT radius R, so that there exists an opti- 
mal radius i? Q pt — 1-46 nm at which f|l 32[) and (|133[) are 



compatible, 
field is 



For R = i? pt, the compensating magnetic 
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FIG. 16. (Color online) Chiral (23,20)-CNT with electric field 
E — lV/nm. (a) Fhe spectrum and (b) the spin expectation 
values at the K/K' point for the |2, k) subband and magnetic 
field B z = 0. For B z = 0.46T, the bands (c) at K' (not 
shown) are gapped, while the spectrum at K (solid lines) has 
the same form as in the armchair case [Eq. 1|126[) ]. The size 
of the gap is 2\a\ = 0.16meV. The spin expectation values (d) 
at K for |2, k) follow closely the armchair case [see Eq. I)129p ]. 



B z>opt = 3.67T-cos30. (134) 

Of course, the CNT radius is not a continuous variable, 
but can only take on discrete values. Fig. [15] shows 
the fields needed for compensating the effective Zeeman 
term for different CNT chiralities. We see in Fig. [T5"lthat 
there are indeed CNTs that have an optimal radius. In 
fact, the uncertainty in determining the optimal radius 
of actual CNTs is larger than the spacing of the dots in 
Fig. 14 so that all dots in the vicinity of the vertical line 
in the figure can be considered as optimal. 
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FIG. 15. (Color online) Magnetic fields needed for compen- 
sating the chirality-induced effective Zeeman term in different 
(m + 31, m)-CNTs as a function of CNT radius R. I varies 
from 1 (lower curves, red) to 5 (upper curves, blue). The dots 
correspond to the discrete radii and chiralities. The intercon- 
necting lines are guides to the eye. The vertical line shows 
the optimal CNT radius at which both, the Zeeman term and 
the orbital shift, are compensated simultaneously. 



D. External magnetic fields in armchair CNTs 

In this section we discuss the low-energy spectra of 
armchair carbon nanotubes in strong electric fields and 
additional magnetic fields B along the three possible spa- 
tial directions. 



1. Magnetic field along the nanotube 

A magnetic field along the CNT axis B = B z z leads 
to the shift of the circumferential wave vector, given by 
Eq. (|13ip . Together with the Zeeman term, this leads to 
the additional term in the Hamiltonian 

#mag = nvp^—ax + fi B B z S z . (135) 

The resulting spectrum of a (20,20)-CNT in a 1 V/nm 
electric field is shown in Fig. [17] for different magnetic 
field strengths. 

As explained above, for the case of a magnetic field 
along the CNT axis, the additional terms in the Hamil- 
tonian can be accounted for by redefining the parameters 
/3 and Afc* v , i.e., 

f3*=f3 + fi B B z , (136) 
Afc** = A*4 + ttB z R/3> . (137) 

The energy spectrum at k = for a chiral nanotube is 
given by 

£1.3 = -a ± v/(e££) 2 + W* - hvpAk**) 2 , (138) 
£ 2 .4 = a± ^J{eE£) 2 + (/3* + hvpAk**) 2 . (139) 
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The splitting Ai at k — is given, to leading order in 
the magnetic field B z and for an armchair nanotube, by 



(/3* +hv F Ak t *Y 
2eEi 



(140) 
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FIG. 18. (Color online) The spectrum of a (20,20)-CNT in the 
electric field E = lV/nm and different magnetic fields B x = 
0,0.5, 1 T (parts a,b,c) along the electric field (x direction). 
The dependence of the spin-polarization on the wave vector k 
is represented on the graph d) for B x = 0.1 T. In the limit of 
large magnetic fields the Zeeman term dominates and (S x ) — > 
±1. See Eq. (fll3]l for Ai. 



FIG. 17. (Color online) Spectrum of a (20,20)-CNT in 
an electric field E = IV /nm and different magnetic fields 
B z = 0, 0.09, 0.03 T (parts a,b,c) along the nanotube axis (z 
direction). The dependence of the spin-polarization on the 
wave vector k is shown in part d) for B~ — 0.03 T . In the 
limit of the large magnetic fields the Zeeman term dominates 
and (S z ) -t ±1. See Eq. (fT40|) for Ai. 



2. Magnetic field along the electric field 

In case of a magnetic field perpendicular to the CNT 
axis we take into account only the Zeeman term, as the 
orbital effect is small for the strengths of the fields con- 
sidered here. For B = B. x ic parallel to the electric field 
the additional term in the Hamiltonian reads 

Hmag = fl B B x S X . (141) 

For an armchair CNT, the total Hamiltonian can then 
be diagonalized analytically. We find 




-0.4 0.4 -0.4 0.4 



k(fim x ) k(/j,m 1 ) 

FIG. 19. (Color online) The spectrum of the CNT (20,20) 
in the electric field E = IV /nm and different magnetic fields 
By — 0, 0.2, 0.5, 1 T (parts a,b,c,d) perpendicular to the elec- 
tric field and to the nanotube axis (y direction). The depen- 
dence of the spin-polarization on the wave vector is given by 
Eq. (|129[) with the shifted wave vector. 



e = ± 



(e£0 2 + « 2 + {hbB x ) 2 + [hv F kf 



± 2 ^(e^a) 2 + ((e££) 2 + (^bB x ) 2 ) (hv F k) 2 

"(142) 

In Fig. [T5]the spectrum is shown for several values of the 
magnetic field. B x decreases the splitting at k = 



Ai ~ 2 



, E f (a 1 \ 
K V 2 a 2 - (e^) 2 J 



(143) 



3. Magnetic field perpendicular to the electric field and the 
nanotube axis 



A magnetic field perpendicular to the electric field and 
to the nanotube axis is described by the Zeeman term 



H mag — ^bB v S v . 



(144) 



In this case the magnetic field tries to polarize the spin in 
the same direction as the electric field. For an armchair 
nanotube the Hamiltonian can be diagonalized exactly 
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giving the energy spectrum 



£1,2 = -e££ ± J a 2 + (Hv F k - (J,bB v ) 2 , (145) 



£3.4 = eE£ ± Jo? + (hvFk + hbB. 1)2 



VI > 



(146) 



which is shown in Fig. [19J B y leads to an additional shift 
Afc = dziigBy/hvp of the momentum along the nanotube 
axis. 



VII. RESONANT SPIN TRANSITIONS 



leads to an additional complication compared to a simple 
spin resonance Hamiltonian of the form cos(ujt)S x . This 
complication disappears in the limit of large fc where the 
term rfrupfca^ dominates in and one can assume 02 
to be a good quantum number, i.e., u<i = ±1. In this 
case the Rabi frequency is given by 



UJr 



eE ac i 



(149) 



The occurence of the sublattice operator in the coupling 
term reduces the Rabi frequency uir < lo r around k ~ 0. 



The manipulation of the electron spin by time- 
dependent external fields is most important for quan- 
tum computing and spintronics. Traditionally, a time- 
dependent magnetic field, coupling to the electron spin 
via the Zeeman energy, is utilized for Rabi flopping^ 
However, the combination of time-dependent electric 
fields and spin-orbit coupling may give rise to an all- 
electric control of the electron spin42r— This effect is 
called electric dipole spin resonance (EDSR) and it was 
argued semiclassically^ that this concept is applicable to 
CNTs in principle. In the following, we investigate the 
EDSR effect due to Hf Ql which has been derived micro- 



Without dc electric field 



scopically in Sec. IIVI 

In a spin resonance experiment one considers a spin 
which is split by a Zeeman energy Ez- In order to drive 
transitions between the spin-up and spin-down states, 
an external field with frequency u = Ez /h is required. 
In the case of a CNT, there are two regimes for such 
transitions, characterized by the frequencies required to 
drive them. Spin resonance at optical frequencies (THz 
regime) involve different subbands. They arc possible^ 
but not of interest in this work. The intra-subband tran- 
sitions have characteristic energy scales below 1 meV, 
which corresponds to frequencies in the GHz regime. AC 
voltages in the GHz regime can easily be generated elec- 
tronically, so that this regime is suitable for EDSR. 

We now discuss the transitions between electronic 
states with opposite spins in the lowest subband, induced 
by a time-dependent electric field 



E ac (t) = EacCOSLOt. 



(147) 



The effective Hamiltonian describing the interaction with 
a time-dependent field is similar to Eq. (|74|) . However, 
we assume E ac along the y direction, so that 



H ac (t) = TeE ac {t)iS x a 2 . 



(148) 



This term can be used to implement the EDSR effect. 
The frequency w of the field is chosen to fit the energy 
difference between two eigenstates of H° s [Eq. (JBSJ], 
between which the transitions are induced. 

Indeed, the form of H ac (t) is typical for the EDSR 
effect. It has a trigonometric time dependence and is 
proportional to the spin operator S x . However, H &c {t) is 
also proportional to the sublattice operator 02, and this 



We start with considering non-armchair CNTs with 
9 7^ 7r/6. In this case, transitions between the states e u + 
and e u - [see Eq. (|12ip ] are driven by H ac {t). For the 
Rabi frequency we obtain 



OJR 



eE ac £ 



sin((t?++i?_)/2)| 



(150) 



where i!}± is defined in Eq. (|122[) . In agreement with 
what was explained above, Eq. (|150[) reduces to u* R in 
the limit k — > ±00. 

For k = 0, the wavefunction is an eigenfunction of 
(j 1 . The states £„,+ and e u - have opposite spins but 
the same isospins a\, if HvF(k t + Afc* v ) > a. The time- 
dependent electric field [Eq. (|148j) ]. however, couples spin 
(S) and isospin (a) simultaneously. Thus, at k = Rabi 
flopping between e„ j+ and e Ut - via H ac (t) is not allowed 
and the Rabi frequency is zero (see Fig. [20)) . Near the 
Dirac points the Rabi frequency is proportional to k. We 
find 



eE ac i 



fc* + Ak f r 



(151) 



For non-metallic CNTs, kt is very large and this leads to 
a strong suppression of the Rabi frequency (see Fig. |2"0"1) . 

A magnetic field along the nanotube axis renormalizes 
the coefficients (3 and Afc' v [see Eqs. (fT36)) and (fT37)) ] 
and thus allows us to change the Rabi and resonance 
frequencies. If the magnetic field is chosen such that 
fe* + Afc** = 0, i.e., for B = Bf [see Eq. (|T^2")) ]. with 



Bt 



<f>o 5.4meV 
irhvF i?[nm] 3 



cos W. 



(152) 



the Rabi frequency is u* R for arbitrary fc. Thus, by 
chosing B z appropriately, the Rabi frequency can be in- 
creased to its upper limit lo* r . This effect is stable with 
respect to small deviations from the optimal B z , as is 
shown in Fig. Note that for fc = 0, lor = uj* r for any 
B z sufficiently close to B z r . 

As usual, a magnetic field B perpendicular to the nan- 
otube axis aligns the electron spin. An ac electric field 
along this magnetic field induces Rabi transitions. As- 
suming that the frequency of the electric field is tuned 
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FIG. 20. (Color online) The dependence of the Rabi frequency ujr on the wave vector k along the CNT axis (a) for the 
metallic nanotubes (13,10)-CNT (dotted), (23,20)-CNT (dashed), and (33,30)-CNT (solid), (b) shows ujr for semiconducting 
nanotubes: the CNT (12,10) - dotted, the CNT (18,10) - dashed, the CNT (24,10) - solid. The amplitude of the ac electric 
field is 1 mV/nm. 
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FIG. 21. (Color online) The dependence of the Rabi frequency 
ujr on the wave vector k along the CNT (26,20) axis. The 
amplitude of the ac field is 1 mV/nm. The magnetic field 
along the CNT axis is equal to B — 0.93Bf (dash-dotted), 
B = 0.95-Bf (dashed), B = 0.9857 (dotted), and B = Bf 
(solid). 
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FIG. 22. (Color online) The dependence of the Rabi frequency 
ujr on the magnetic field B± applied perpendicular to the 
CNT axis for metalic nanotubes: the CNT (13,10) -dotted, 
the CNT (23,20) - dashed, the CNT (33,30) - full. The am- 
plitude of the ac field is 1 mV/nm. 



With static electric field 



to the energy of the spin splitting at the Fermi points 
(the Fermi level is assumed to be tuned into the k = 
splitting), the dependence of the Rabi frequency on the 
magnetic field is shown in Fig. 

If a circularly polarized field is applied, one induces the 
transitions only from the spin-down state at the Fermi 
level to the state spin-up above the Fermi level 

H circ (t) = TeE clrc (t)£S+a 2 , (153) 

with = S x ± iS y . If the opposite polarization of 
the field is chosen, the transitions occur in the opposite 
directions. 



A static electric field perpendicular to the axis of a nan- 
otube in combination with SOI lifts the spin degeneracy 
in the spectrum of an armchair nanotube. In contrast 
to the well-known Rabi resonance method, realized in a 
static magnetic field and a perpendicular time-dependent 
magnetic field, we propose an all-electric setup for spin 
manipulation with two perpendicular electric fields, one 
of which is static and the other is time-dependent. The 
static electric field aligns the spin along the y direction, 
i.e., perpendicular to the CNT axis and perpendicular 
to the direction of the static electric field. The time- 
dependent electric field rotates spin around the x direc- 
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tion. 

For armchair nanotubes we find that the transitions 
between states 1 and 2 or between 3 and 4 in the spec- 
trum shown in Fig. UOI have the optimal Rabi frequencies 
wjj. Rabi transitions between these groups (e.g. 1h4) 
are not possible. 
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FIG. 23. (Color online) The dependence of the Rabi frequency 
ujr on the wave vector k along the CNT (26,20) axis. The 
amplitude of the ac field is E ac = lmV/nm and E — IV /nm. 
The magnetic field along the CNT axis is equal to B = 
(dash-dotted), B = 0.5B C / (dashed), B = 0.9B C / (dotted), 
and B = B% r (solid). 

As explained above, for non-armchair but metallic 
CNTs, the orbital term Afc** can be compensated by 
an additional magnetic field B z = B c z r along the CNT. 
Then, as in the case of the armchair nanotubc, the Rabi 
frequency of transitions between states I and 2 or be- 
tween 3 and 4 is uj* r , while transitions between states of 
different groups are not allowed. For B z ^ B z r , the Rabi 
frequencies are smaller than lo* b and depend on k, as is 
shown in Fig. [23l 

A magnetic field perpendicular to the nanotube axis 
(see FigfT8j) breaks the symmetry of the spectrum around 
the Dirac point. As a result, the resonance frequencies 
for the right-moving and left-moving modes are different 



and it is possible to implement the EDSR mechanism for 
only one of the two modes. 



VIII. CONCLUSIONS 

We have studied the interplay of strong electric fields, 
magnetic fields and spin-orbit interactions in carbon nan- 
otubes. An approximate effective low-energy theory de- 
scribing the electrons near the two Dirac points has 
been derived analytically and this theory has been tested 
against more sophisticated numerical solutions of the lat- 
tice tight-binding Hamiltonian for the second shell 7r and 
a orbitals. We have established that the properties of 
carbon nanotubes are described well by our analytical 
model in the limit of large field screening. The latter can 
be achieved by immersing the CNT into dielectrica. 

The central feature of CNTs in electric fields is the ap- 
pearance of (spin-filtered) helical modes in an all-electric 
setup. For perfect armchair nanotubes, there are two 
pairs of helical modes, one for each valley, transporting 
up-spins in one direction and down-spins in the opposite 
direction. This helicity is perfect in that the average spin 
is non-zero only for this one spin component and zero for 
all others. Thus, the average spin is a perfectly odd func- 
tion of k for armchair CNTs. For non-armchair chirali- 
ties, an additional magnetic field can be used to restore 
the helical phase in one valley. In the other valley, all 
electronic states are removed from the low-energy regime 
so that this valley is suppressed by the combination of 
non-armchair chirality and magnetic field. 

Furthermore, we have shown that the EDSR effect 
may be implemented by a time-dependent electric field 
perpendicular to the CNT. The typical Rabi frequencies 
which can be achieved in this system are in the MHz-GHz 
range. 
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